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Abstract 

In this paper, some properties of the minimal speeds of pulsating Fisher-KPP fronts 
in periodic environments are established. The limit of the speeds at the homogenization 
limit is proved rigorously. Near this limit, generically, the fronts move faster when the 
spatial period is enlarged, but the speeds vary only at the second order. The dependence 
of the speeds on habitat fragmentation is also analyzed in the case of the patch model. 



1 Introduction and main hypotheses 

In homogeneous environments, the probably most used population dynamics reaction- 
diffusion model is the Fisher-KPP model [13, 23] . In a one-dimensional space, it corresponds 
to the following equation 

— ^D— + u (i^-iyu), t>Q, X eR. (1.1) 

The unknown u — u{t, x) is the population density at time t and position x, and the positive 

constant coefficients D, n and u respectively correspond to the diffusivity (mobility of the 
individuals), the intrinsic growth rate and the susceptibility to crowding effects. 

A natural extension of this model to heterogeneous environments is the Shigesada- 
Kawasaki-Teramoto model [32], 

Ozt I Ou \ 

di^dx V^^^^dx) il^Lix) -lyMu), t > 0, x eW, (1.2) 



where the coefficients depend on the space variable a; in a L-periodic fashion: 



Definition 1.1 (L-periodicity) Let L be a positive real number. We say that a function 
/i : R — > R is L-periodic if 

V x e R, h{x + L) = h{x). 
In this paper, we are concerned with the general equation: 



du d f , ,du 
dt dx 



aL(x) — j +fL{x,u), i e R, X e R. (1.3) 



The diffusion term aL satisfies 

ciL{,x) = a{x/L), 

where a is a C^'^{R) (with S > 0) l-periodic function that satisfies 

3 < ai < q;2, V x e R, ai < a{x) < a^. (1.4) 

On other hand, the reaction term satisfies //.(a;, ■) = f (x/L,-), where f := f{x,s) : RxR+ — >• 
R is l-periodic in x, of class C^'^ in {x, s) and in s. In this setting, both and are 
L-periodic in the variable x. Furthermore, we assume that: 

VxgR, /(a;,0)=0, 

3M>0, Vs>M, VxeR, fix,s)<Q, (1.5) 
V X e R, s I— > f{x, s)/s is decreasing in s > 0. 



Moreover, we set 
and 



H{x) := lim f{x,s)/s, 

s— *0+ 



:= ^lim /^(x, s)/s = /i (^|) 



The growth rate // may be positive in some regions (favorable regions) or negative in others 
(unfavorable regions). 

The stationary states p{x) of (1.3) satisfy the equation 

^ («^(^)^) + /^(^'^') = 0, ^ e R. (1.6) 

Under general hypotheses including those of this paper, and in any space dimension, it was 
proved in [4] that a necessary and sufficient condition for the existence of a positive and 
bounded solution p of (1.6) was the negativity of the principal eigenvalue pi^L of the linear 
operator 

Co: ^^-iaL{xWy-l^L{x)^, (1-7) 

with periodicity conditions. In this case, the solution p was also proved to be unique, and 
therefore L-pcriodic. Actually, it is easy to sec that the map L pix is nonincrcasing 
in L > 0, and even decreasing as soon as a is not constant (see the proof of Lemma 3.1). 
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Furthermore, pi^i ^ ~ l^{x)dx as L ^ 0^. In this paper, in addition to the above- 

Jo 

mentioned hypotheses, we make the assumption that 



1 

li{x)dx > 0. (li 



'0 

This assumption then guarantees that 

V L > 0, < 0, 

whence, for all L > 0, there exists a unique positive periodic and bounded solution pl of 
(1.6). Notice that assumption (1.8) is immediately fulfilled if n{x) is positive everywhere. 

In this work, we are concerned with the propagation of pulsating traveling fronts which 
are particular solutions of the reaction-diffusion equation (1.3). Before going further on, we 
recall the definition of such solutions: 

Definition 1.2 (Pulsating traveling fronts) A function u = u(t,x) is called a pulsating 
traveling front propagating from right to left with an effective speed c ^ 0, if u is a classical 
solution of: 

VA; e Z, V(t,a;) e M X M, u(t^ ,x) ^ u(t,x + kL), (1.9) 

< u{t,x) < Pl{x), 
lim u{t,x) = and lim u{t,x) — p]:,{x) = 0, 

x—*—oo x—*+oo 

where the above limits hold locally in t. 

This definition has been introduced in [31, 32]. It has also been extended in higher 
dimensions with pi, = 1 in [1] and [35], and with ^ 1 in [5]. 

Under the above assumptions, it follows from [5] that there exists > such that 
pulsating traveling fronts satisfying (1.9) with a speed of propagation c exist if and only 
if c > c^. Moreover, the pulsating fronts (with speeds c > c^) are increasing in time t. 
Further uniqueness and qualitative properties are proved in [14, 15]. The value is called 
the minimal speed of propagation. We refer to [2, 3, 11, 18, 25, 27, 28, 34] for further 
existence results and properties of the minimal speeds of KPP pulsating fronts. For existence, 
uniqueness, stability and further qualitative results for combustion or bistable nonlinearities 
in the periodic framework, we refer to [6, 7, 12, 16, 17, 19, 24, 26, 35, 36, 37, 38]. 



In the particular case of the Shigesada et al model (1.2), when a{x) = 1, the effects of 
the spatial distribution of the function /^l on the existence and global stability of a positive 
stationary state pl of equation (1.2) have been investigated both numerically [30, 31] and 
theoretically [4, 8, 29]. In particular, as already noticed, enlarging the scale of fragmentation, 
i.e. increasing L, was proved to decrease the value of pi^^. Biologically, this result means 
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that larger scales have a positive effect on species persistence, for species whose dynamics is 
modelled by the Shigesada et al model. 

The effects of the spatial distribution of the functions ai and fiL on the minimal speed of 
propagation have not yet been investigated rigorously This is a difficult problem, since 
the known variational formula for bears on non-self- adjoint operators, and therefore, 
the methods used to analyze the dependence of pi^L on fragmentation cannot be used in this 
situation. However, in the case of model (1.2), when gl = 1, = 1 and [Il^x) = fi{x/L), for 
a 1-periodic function fi taking only two values, Kinezaki et al [22] numerically observed that 

was an increasing function of the parameter L. For sinusoidally varying coefficients, the 
relationships between and L have also been investigated formally by Kinezaki, Kawasaki, 
Shigesada [21]. The case of a rapidly oscillating coefficient ai(x), corresponding to small 
L values, and the homogenization limit L ^ 0, have been discussed in [19] and [38] for 
combustion and bistable nonlinearities f{u). 

The ffist aim of our work is to analyze rigorously the dependence of the speed of propaga- 
tion with respect to L, under the general setting of equation (1.3), for small L values. We 
determine the limit of the minimal speeds as L — > 0"*" (the homogenization limit) , and we 
also prove that near the homogenization limit, the species tends to propagate faster when the 
spatial period of the environment is enlarged. Next, in the case of an environment composed 
of patches of "habitat" and "non-habitat" , we consider the dependence of the minimal speed 
with respect to habitat fragmentation. We prove that fragmentation decreases the minimal 
speed. 



2 Main results 

In this section, we describe the main results of this paper. Unless otherwise mentioned, we 
make the assumptions of Section 1. The first theorem gives the limit of as L goes to 0. 

Theorem 2.1 Let c*j^ be the minimal speed of propagation of pulsating traveling fronts sol- 
ving (1.9). Then, 

lim cl = 2^<a>H <H>a, (2.1) 

where 



<^>A — j iJ>{x)dx and <a>H — {^j {(^ix)) ^dx^ 



denote the arithmetic mean of n and the harmonic mean of a over the interval [0, 1]. 

Formula (2.1) was derived formally in [33] for sinusoidally varying coefficients. Theo- 
rem 2.1 then provides a generalization of the formula in [33] and a rigorous analysis of the 
homogenization limit for general difi^usion and growth rate profiles. 

Remark 2.2 The previous theorem gives the limit of as L — when the space dimension 
is 1. Theorem 3.3 of El Smaily [11] answered this issue in any dimensions N, but under an 
additional assumption of free divergence of the diffusion field (in the one-dimensional case 
considered here, this assumption reduces to da/dx = in M). Lastly, we refer to [6, 7, 16] 
for other homogenization limits with combustion-type nonlinearities. 
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Our second result describes the behavior of the function L i— > c^, for small L values. 



Theorem 2.3 Let c*j^ be the minimal speed of propagation of pulsating traveling fronts sol- 
ving (1.9). Then, the map L ^ is of class C°° in an interval (0, Lq) for some Lq > 0. 
Furthermore, 

rlr* 

lim ^ = (2.2) 

and 

lim = 7 > 0. (2.3) 

L^o+ dL^ 



Lastly, J > if and only if the function 



fJ- _^ <a>H 



<IJi>A 

is not identically equal to 2. 

CoroUciry 2.4 Under the notations of Theorem 2.3, it follows that if a is constant and /i 
is not constant, or if /i is constant and a is not constant, then 7 > and the speeds are 
increasing with respect to L when L is close to 0. 

Remark 2.5 The question of the monotonicity of the map L 1-^ had also been studied 
under different assumptions in [11] (see Theorem 5.3). The author answered this question for 
a reaction-advection-diffusion equation over a periodic domain Q C M^, under an additional 
assumption on the diffusion coefficient (like in Remark 2.2, this assumption would mean 
again in our present setting that the diffusion coefficient a{x) is constant over R). Our result 
gives the behavior of the minimal speeds of propagation near the homogcnization limit for 
general diffusion and growth rate coefficients. The condition 7 > is generically fulfilled, 
which means that, roughly speaking, the more oscillating the medium is, the slower the 
species moves. But the speeds vary only at the second order with respect to the period L. 
Based on numerical observations which have been carried out in [21] for special types of 
diffusion and growth rate coefficients, we conjecture that the monotonicity of holds for 
all L > 0. 

Lastly, we give a first theoretical evidence that habitat fragmentation, without changing 
the scale L, can decrease the minimal speed c*. We here fix a period Lq > 0. 

We assume that a = 1, and that /i^g := /i^ takes only the two values and m > 0, and 
depends on a parameter z. More precisely: 

There exist < z and / G (0, Lq) such that I + z < Lq, 
IJ^=mon[0,l/2)U[l/2 + z,l + z), (2.4) 
/i, = on [1/2, 1/2 + z)U[l + z, Lq). 

With this setting, the region where ji^ is positive, which can be interpreted as "habitat" in 
the Shigesada et al model, is of Lebesgue measure / in each period cell [0, Lq]. For z = this 
region is simply an interval. However, whenever z is positive, this region is fragmented into 
two parts of same length 1/2 (see Figure 1). Our next result means that this fragmentation 
into two parts reduces the speed c*. 
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(a) 



l+z 



(b) 



Figure 1: The Lo-periodic function x ^— /U^(a;), (a): with 2; = 0; (b): with z > Q. 



Theorem 2.6 Let c* he the minimal speed of propagation of pulsating traveling fronts sol- 
ving (1.9), with a^^ = 1 and = /i^ defined by (2.4). Assume that I G (3Lo/4, Lq). T/ien 
2; I— i> c* decreasing in [0, (Lq — /)/2], anc? increasing in [{Lq — l)/2, Lq — I]. 

Remark 2.7 Note that, whenever z > (Lq — /)/2, the two habitat components in the period 
cell [1/2 + z,Lq + 1/2 + z] are at a distance smaller than (Lg — l)/2 from each other. In 
fact, Theorem 2.6 proves that, when z varies in (0, L^ — l), c* is all the larger as the minimal 
distance separating two habitat components is small, that is as the maximal distance between 
two consecutive habitat components is large. 

Remark 2.8 Here, the function fiz does not satisfy the general regularity assumptions of 
Section 1. However, c* can still be interpreted as the minimal speed of propagation of weak 
solutions of (1.9), whose existence can be obtained by approaching fi^ with regular functions. 

The main tool of this paper is a variational formulation for involving elliptic eigenvalue 
problems which depend strongly on the coefficients a and /. Such a formulation was given in 
any space dimension in [3] in the case where the bounded stationary state p of the equation 
(1.3) is constant, and in [5] in the case of a general nonconstant bounded stationary state p{x). 



3 The homogenization limit: proof of Theorem 2.1 

This proof is divided into three main steps. 

Step 1: a rough upper bound for c^. For each L > 0, the minimal speed is positive 
and, from [5] (see also [3] in the case when p = 1), it is given by the variational formula 

^ A>o A A* ' ^ ' 

where A^ > and, for each A G M and L > 0, A;(A, L) denotes the principal eigenvalue of the 
problem 

(«lV^a,l)' + '^^(^Li^'x,L + ^olhi^KL + A^olV^a.l + /WlV'a,l = k{\ L)i^\L in M, (3.2) 

with L-periodicity conditions. In (3.2), ipx^L denotes a principal eigenfunction, which is of 
class C^''^(]R), positive, unique up to multiplication by a positive constant, and L-periodic. 
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Furthermore, it follows from Section 3 of [5] that the map A i— > A;(A, L) is convex and that 
II (0, L) = for each L > 0. Therefore, for each L > 0, the map A i— > /c(A, L) is nondecreasing 
in M_|_ and 

V A > 0, V L > 0, /c(A, L) > k{0, L) = > (3.3) 

under the notations of Section 1. 

Multiplying (3.2) by ipx^^ and integrating by parts over [0,L], we get, due to the L- 
periodicity of and ■?/'a,l ■ 

rL i*L pL pL 

Jo Jo Jo Jo 

for all A > and for all L > 0. Consequently, 

VA>0, VL>0, k{X,L) < X'^ttM + ^iM, (3.4) 

where 

qm — maxa(x) > and um — max/i(x) > 0. 
Using (3.1), we get that 

V L > 0, 0<cl< 2,/^]^. (3.5) 
Step 2: the sharp upper bound for c^. For any A > and L > 0, consider the functions 

Since ipx^L is unique up to multiplication, we will assume in this step 2 that 

!\lAx)dx^l. (3.6) 

The above choice ensures that 

/ i,lL{x)dx< [ e''-^2PlL{x)dx= [ ^lr.{^)dx = l. (3.7) 
JO Jo Jo 

We are now going to prove that the families {iPx,l)\,l and (93a,l)a,l remain bounded in 
H^{0, 1) for L small enough and as soon as A stays bounded. For each L > 0, we call 

Ml = [l/L] + 1 e N, 

where [l/L] stands for the integer part of l/L. Multiplying (3.2) by V'a.l and integrating by 
parts over [0, Ml L], we get that 



pM^L pMlL i-MlL i-MlL 

/ aLi^'x,L+ AW'A,L+/ f^L^lL-k{X,L) -2 

Jo Jo Jo Jo 



X,L- 
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Using (1.4), (3.3) and (3.4), it follows that 

< / 1/^'^/ < - X (A'oM + //m) X / Va,l- 

Jo 0(1 Jo 

Since 1 < MlL < 1 + L for all L > 0, we have that 1 < MlL < 2 for all L < 1. Thus, for 
all < L < 1, 

/ i^'x,L < / ^'x/ and / V'U ^ / < 1 

Jo Jo Jo Jo 

from (3.7). It follows now that 



VA>0, V0<L<1, /'Vl/< ^'"^ + ^^ . (3.8) 

Jo ' «i 

Prom (3.7) and (3.8), we conclude that, for any given A > 0, the family (V'a,l)o<a<a, o<l<i 
is bounded in H^{0, 1). On the other hand, 

^'x,l{^) = Av5A,L(a;) + e^^^l,i(a;). 
Owing to (3.6) and (3.8), we get: 

VA>0, VL<1, \\^'x,l\\lho,i) < a ||v2a,l||l2(o,i) +e^||V'l,LllL2(o,i) 



(3.9) 



< A + X 



A^om + Mm 



Oil 



From (3.6) and (3.9), we obtain that, for any given A > 0, the family ((^a,l)o<a<a, o<l<i is 
bounded in H^{0,1) and that the family {aL^'^^ l)o<x<a, o<l<i is bounded in L^(0, 1) (due 
to (1.4)). Moreover, 

(aL<^A,L)' = A^az.e^^V'A.L + 2Aai.e^-V'l + Aa'^e^^V^^.L + e^'a'^V'l.L + e'^a^.^^,^. 
Multiplying (3.2) by e'^^, we then get 

{(^L^'x,Ly + I^L^x,L = k{\ L)(px,L in R. (3.10) 

Let 

for all A > 0, L > and a; e M. Pick any A > 0. One already knows that the family 
(^a,l)o<a<a, o<l<i is bounded in 1/^(0, 1). Furthermore, 

^A,L + t^L<Px,L = k{\ L)cpx,L in M. (3.11) 

Notice that the family (A;(A, L))o<a<a, o<l<i is bounded from (3.3) and (3.4). Prom (3.6) 
and (3.11), it follows that the family (^a,l)o<a<a o<l<i bounded in -L^(0, 1). Eventually, 
(^^a,l)o<a<a, o<l<i is bounded in H^{0, 1). 
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Pick now any sequence (L„)„gj^ such that < L„ < 1 for all n e N, and L„ — > 0+ as 
n — i> +00. Choose any A > and any sequence (A„)„gN of positive numbers such that A„ — > A 
as n —>■ +00. We claim that 

k(Xn,Ln)^X'^<a>H + <lJi'>A as n ^ +00, (3-12) 
where <a>H— {^j {a{x))~^dx^ and <^>a= j n{x)dx. To do so, call 

■Ipn = ^A„,L„, (fn = <PX„,L„ and = Vx„,L„. 

It follows from the above computations that the sequences and (vn) are bounded in 
H^{0,1). Hence, up to extraction of a subsequence, 

ipn ^ and — > w as n — > +oo, 

strongly in L^(0, 1) and weakly in H^{0,1). By Sobolev injections, the sequence (ipn) is 
bounded in C°'^/^([0, 1]). But since each function ■0„ is L„-periodic (with L„ — > 0+), it follows 
from Arzela-Ascoli theorem that ijj has to be constant over [0, 1] . Moreover, the boundedness 
of the sequence {k{Xn, i^n))„gN ™plies that, up to extraction of another subsequence, 

/c(A„, Ln) — >■ k{X) e R as n ^ +oo. 

We denote this limit by A; (A), we will see later that indeed it depends only on A. It follows 
now, from (3.11) after replacing (A, L) by (A„, L„) and passing to the limit as n — > +oo, that 

w'+ <n>A e^'^if = A-(A)V'e^^ a.e. in (0, 1). 
Notice indeed that /il ^<pi'>A as L — > 0+ in -L^(0, 1) weakly. Meanwhile, 

(p'^ = A^e^^^V^n + e^^^Vn = — ^ <a~^>A w as n ^ +oo, weakly in L^(0, 1), 



where <a >a 



/ (a(x)) ^dx. Thus, we obtain 
Jo 



Consequently, 

A^ < a >ij V+ < A* >A V' = k{X)ip. 
Actually, since the functions ■0„ are L„-periodic (with L„ — > 0^) and converge to the constant 
■0 strongly in L^(0, 1), they converge to -0 in Lfoci^)- 

Jo Jo 

where M = sup„gpj A^. Hence, ip ^ and 

X'^ <a>H + <l^>A^k{X). (3.13) 
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By uniqueness of the limit, one deduces that the whole sequence (/c(A„, Ln))n&n converges to 
this quantity A;(A) as n ^ +00, which proves the claim (3.12). 

Now, take any sequence L„ 0"*" such that c^^ limsupi_^o+ c}^ SiS n —>■ +00. For each 
A > and for each n eN, one has 



* ./c(A,L„) 

Cr < ^ 



from (3.1), whence 



r * 1. * ^ H^) X / ^ , </'>A 
limsupcj^= lim Cj^ < — r — = A <a>H -\ ^ • 

L—^0+ n-*+oo "A A 

Since this holds for all A > 0, one concludes that 

lim sup < 2y/<a>H <I^>a- (3-14) 

Step 3: the sharp lower bound for c^. The aim of this step is to prove that 

lim inf c*r > 2J<a>H <IJ.>a 

which would complete the proof of Theorem 2.1. 

For each L > 0, the minimal speed is given by (3.1) and the map (0, +00) 9 A 1— > 
A:(A, L)/X attains its minimum at A^ > O.Wc will prove that, for L small enough, the family 
(A^) is bounded from above and from below by A > and A > respectively. Namely, one 
has 

Lemma 3.1 There exist Lq and < A < A < +00 such that 

X < XI <X for all < L < Lq. 

The proof is postponed at the end of this section. Take now any sequence {Ln)n such 
that < Ln < Lq for all n, and L„ ^ 0^ as n ^ +00. From Lemma 3.1, there exists A* > 
such that, up to extraction of a subsequence, A^^ — > A* as n — > +00. One also has 

k{K„,Ln) k{X*) ^ , </^>A ^ ^ . ^ ^ ^ 

= — ^ = A <a>H > 2V<a>H<f^>A 

J-fn 

from (3.12) and (3.13). Therefore, hminfi_^o+c^ > '^\/<cl>h <pi'>A- Eventually, 



lim^ = 2^<a>H <I^>a 

and the proof of Theorem 2.1 is complete. □ 

Proof of Lemma 3.1. Observe first that, for A = and for any L > 0, k{0,L) is the 
principal eigenvalue of the problem 

{aM' + I^L(I>L = k{0, L)(t)L in K, 
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and we denote 0l = ipo,L a principal eigenfunction, which is L-periodic, positive and unique 
up to multipUcation. In other words, k{0, L) = — pi,L under the notations of Section 1. 
Dividing the above eUiptic equation by 0l and integrating by parts over [0, L], one gets 

1 /-^ai^'/ 



k{0,L) = - —7^+ i^{x)dx > <i^>A > 0. 

^ Jo 0L Jo 

On the other hand, as already recalled, ||(0, L) — and the map A i— > k{X, L) is convex for 
all L > 0. Therefore, 

V A > 0, V L > 0, k{X,L) > k{0,L) > <h>a > 0. 

Assume here that there exists a sequence (L„)„gN of positive numbers such that L„ — > 0"*" 
and Ar ^ 0+ as n — > +oo. One then gets 

Cl„ = Y* ~T* +00. 

This is contradiction with (3.14). Thus, for L > small enough, the family {X*i)l is bounded 
from below by a positive constant A > (actually, these arguments show that the whole 
family {\i)l>q is bounded from below by a positive constant). 

It remains now to prove that is bounded from above when L is small enough. 

We assume, to the contrary, that there exists a sequence L„ — > 0"*" as n — > +oo such that 
Ar —>■ +00 as n — > +oo. Call 

kn = ^(A^„,-L„), ^|Jn{x) = il^xi^,Ln(x) and (^„(x) = iPxi^,lAx) = e^*Ln'''ij;^(x) 

for all n e N and x e R. Rewriting (3.10) for A = A^^ and for L — Ln, one consequently 
gets 

V n e N, {aL,,v'J + tiLn<Pn = KVn in K. (3.15) 
Owing to the positivity and the L„-periodicity of the C^(R) eigenfunction it follows that 

V n e N, 3 e [0, L„], Vn(6'n) = maxt/'„(a;) = max ^„(a;), 

xeK xe[o,Ln\ 

whence 

VnGN, V;(^n) = 0. 
For each n e N, let M^^ = [1/L„] + 1 e N. Thus, 

V n e N, if'^iOn + Mi„L„) = Al„e^^n(^"+^^n^")v.„(^„). 

Multiplying (3.15) by ipn and integrating by parts over the interval [OmOn + Ml^L„], one 
then obtains 

aL^On + L„)(^;(^„ + Mi„ K)^n{On + L„) - aL„(^„)(^;(^„)(^„(^„) 

" V ' 

A(n) 



— s/ 

B{n) C(n) 
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But, for each n e N, Ml^ e N while a^^ and ipn a-re I/„-periodic. Hence, aL^{9n + M^^ L„) = 
aiM: M^n + Ml„ Ln) = V'n(^n), and V'U^n + L„) = V;(^n) = 0. Then, 

> — xA2„V'^(^„)e2^i«(^"+^^"-^") (ai>Oisgivenby (1.4)), 

whenever n is large enough so that 2 < e^^^-n^^"^" (remember that A^^ — > +00 as n — > +00, 
by assumption). Meanwhile, for all n e N, 



\C{n)\ < 



e^^*^n-,pl{x)dx < /.oo X X e2^in(«"+^^n^"), (3.18) 

ZX T 



'On 

where /Xoo = Hiax^^gK |/i(a;)|. On the other hand, (3.1) and (3.5) yield 
for all n e N, whence 

pdn + MLn Ln rdn + MLn Ln 

J On ^ On 

< V^H^ X ^niOn) X e'^in(OMnLn) 

Now, the term B{ri) can be estimated as follows 



(3.19) 



— - Pf>n + {J + ^)I^n 
j=0 '^On+jLn 

'^Ln-'^ l-en + (3 + ^)Ln „ 

< ^ a2 e'^^n(^"+0'+i)-^") / «(a;) + A^^^/^n(a;)) rfx [from (1.4)] 

MLn-^ Ln 

= V a2 e2^inM^-+i)M / (V,;(a;) + AI^^„(a;))'rfa; 



j=0 

since is I/„-periodic. One has 

"^n _ pLn /„;./ \ 2 



/ 



«(^)+ALV'n(^))'rf^ <m)^"(|^ + AL) dx. 



We refer now to equation (3.2). Taking X = X*^ , dividing this equation (3.2) by the 
periodic function ipn and then integrating by parts over the interval [0, L„] , we get 

for all neN. Thus, 

^ aL„ + Xl^^ + iiLn^ knL^ < 2^aMl^M x A^^ L„. 



>0 
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Owing to (1.4), it follows that 



Jo V^n / "1 

Putting the above result into B{n), we obtain, for all n e N, 

B{n) < '^^^'^ xAI„L„m) E e^^in(^»+(^-+^)-") 

i=o 

<(5x V.^(^„)e2^in(^''+^^.^«) X (A1„L„ + 1), 
where /3 = i^OL^^Jauii-M jot-y) x C and C is a positive constant such that 

Vx>0, <Cx(x + l). 

Lastly, let us rewrite equation (3.16) as 

VneN, A{n) + C{n)~kn ipl = B{n). 

j$„ 

Together with (3.17), (3.18), (3.19) and (3.20), one concludes that there exists no G N such 
that for n >no, 

- X V'I(^„)e2^in(^"+^^n Lu) (3.21) 

< Px ^l{9n)e'^l^'-^^^-'^-^ X (AI„L„ + 1). 
Divide (3.21) by A2„V'^(^n)e^^^"^^"+^^"-^"^. Then 

y n>no, — - — - < /? X L„ + 



2 2(A^ )2 A^ - V A 



Passing to the limit as n — > +oo, one has L„ — > 0"*" and A^^ — > +oo, whence ai < 0, which 
is impossible. 

Therefore the assumption that A^^ — > +oo as L„ — > 0^ is false and consequently the 
family is bounded from above by some positive A > whenever L is small (i.e. 

< L < Lq). This completes the proof of Lemma 3.1. □ 
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Remcirk 3.2 Prom Theorem 2.1, one concludes that the map (0, +00) 9 L 1— > can be 
extended by continuity to the right at L = 0"*". Furthermore, for any sequence of 
positive numbers such that L„ ^ 0^ as n — >^ +cxd, one claims that the positive numbers A^^ 

given in (3.1) converge to y^<a^^^X/I>A = \J <a~^ >a</^>a as n — > +00. Indeed 

VneN, cl ^^-^^ 

and Lemma 3.1 implies that, up to extraction of a subsequence, A^ — A* > 0. Passing to 
the limit as n — > +00 in the above equation and due (3.13) together with Step 2 of the proof 
of Theorem 2.1, one gets 

2^<a>H<fi>A = = A* <a>H+-j^, 

whence A* = ^ <a>^<ii>A- Since the limit does not depend on any subsequence, one 
concludes that the hmit of A^, as L ^ 0+, exits and 

^hm^ A2 = \J <a>ff^<ijL>A — y/ <a~^ >a< A* >a- 



The shcirp lower bound of lim inf i_»o+ Cl from the homogenized equation. In the 

following, we are going to derive the homogenized equation of (1.3), which will lead to the 
sharp lower bound of liminfj^^o+ c^. However, to furnish this goal we will only consider for 
the sake of simplicity a particular type of nonlinearities among those satisfying (1.5). In fact, 
the following ideas can be generalized to a wider family of nonlinearities which satisfy (1.5), 
but the proof requires technical extra-arguments which will be the purpose of a forthcoming 
paper. 

Por each L > 0, let m^. be a pulsating travelling front with minimal speed for the 
reaction-diffusion equation 

' y{t,x) eRxR, <UL{t + ^,x) ^UL{t,x + L) <1, (3-22) 
lim UL{t,x) = and lim uiit^x) = 1, 

> x— »— 00 a;— »+oo 

where aiix) — a{x/L), a is a C^''^(R) 1-periodic function satisfying (1.4), ji is a C^''^(R) 
positive 1-periodic function and (7 is a C2(M+) function such that y(0) = ^(1) = and 
u I— > g{u)/u is decreasing in (0, -l-oo). Up to a shift in time, one can assume that 

VL > 0, // UL{t,x)dtdx =\. (3.23) 

Ji(o,i)x(o,i) 2 
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For each L > 0, set fL{x,u) :— f{x/L,u) — fj,{x/L)g{u). In this setting, there holds pl = 1- 
Prom standard parabolic estimates, each function ul is (at least) of class C^(R x R). Denote 

VL\t,x) = aL[x)— — [t,x) and WL[t,x) — -^{t,x) m M x M. 

As already underlined, it follows from [1] that wl = ^ > in R x R for each L > 0. 
Under the notations of the beginning of this section, it follows from (1.4) and (3.2) that 
k{\, L) > A^ai + for all L > and A G R, where /im — mm^/i > 0. Hence, > 2^q;i//^ 
for each L > and liminfi^o+ > '^y/cnJI^ > 0. 

We shall now establish some estimates for the functions u^, vl and wl which are inde- 
pendent of L, in order to pass to the hmit as L — > 0"^. Notice first that standard parabohc 
estimates and the (i, x)-periodicity satisfied by the functions ul imply that, for each L > 0, 
oo.x) = and ul{+oo,x) = 1 in Cf„^(R), and ^^^(±00,3;) = in C/q^(R). 

Let k G N\{0} be given. Integrating the first equation of (3.22) by parts over R x 
{—kL,kL), one obtains 

VL>0, [[ f{j,UL)dtdx^2kL. (3.24) 

Multiplying the first equation of (3.22) by ul and integrating by parts over R x {—kL, kL), 
one then gets 

yL>0,kL = -[[ aL{x)(^\ dtdx+ [[ f{j,UL)uLdtdx. {3.25} 

J jRx{-kL,kL) \ OX J J jRx{-kL,kL) ^ 

Notice that the last integral in (3.25) converges because of (3.24) and < f{x/L,UL)uL < 
f{x/L, Ul)- Together with (1.4), one concludes that for each L > 0, the first integral in (3.25) 
converges and 

VL>0, // f?^Vdtdx<^^ 



x{—kL,kL) 
hiL 

at 



dx J CKi 

Multiply the first equation of (3.22) by ^ and integrate by parts over R x {—kL, kL). Since 




one obtains that 

VL>0, // dtdx^ I F{j,l)dx ^2kLx I I g, (3.26) 

J J'Rx{-kL,kL) V / J{-kL,kL) ^ Jo Jo 

where F{y,s) = f{y,T)dT. It follows from the above estimates that for each compact 
subset K of R, 



VO < L < 1, 



xK 



+ ( ^ 

dt J \ dx 



dtdx<C{K), (3.27) 
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where C{K) is a positive constant depending only on K. 

In particular, for each compact 7^ of M and for each L > 0, | I-u^lI |L2(Rxi^) 
Now, differentiate the first equation of (3.22) with respect to t (actually, from the regularity 
of /, the function wl is of class with respect to x). There holds 

dwr, d ( , ^dwi\ ,x. ,, , ^ ^ 

o-lKx)— — + fJ,{Y)g {ul)wl m R X R. 



dt dx \ dx J L 

Multiply the above equation by wl and integrate by parts over R x {—kL, kL). Prom (1.4) 
and (3.26), it follows that 

atax < 



lx{-kL,kL) \ dx J ai 

where rj is the positive constant defined by 

?7 = max/i(a;) max l^fYw)! max |F(x, 1)1 > — — // n( — )g'(uL)wj^dtdx>0. 

xm ue[o,i] xm ^kL J JT^^(^_kL kL) L 

Then, for each compact X C R, there exists a constant C'{K) > depending only on K 
such that ^ 

V0<L<1, // dtdx<C'{K). (3.28) 



xK 



\ dx ) 



Let {Ln)n& be a sequence of real numbers in (0, 1) such that L„ ^ and c|^^ — >■ 
liminfi^o+ > as n — +oo. It follows from (3.27) and the bounds < ml„ < 1 that 
there exists uq in Hl^J^ x M) such that, up to extraction of a subsequence, ul^ —>■ uq strongly 
in Lf^ci^ ^ almost everywhere in R x R, and 

weakly in L]g^{M. x R) as n — > +oo. 







( duo 


duo\ 


[ dt ^ 


dx J 


[ dt ' 


dx J 



Remember that vl„ — aL„^^^^ and < ai < aL„ < «2 for each n e N. Thus, (3.27) 
yields that for each compact X of R and for each n e N, ||t'L„||L2(Rxi<') ^ a2C{K). Further- 
more, (3.22) implies that 

while < f{x/Ln, UL„it, x)) < k in R x R where k — maxR x max[o,i] > is independent 
of n. Together with (3.27), one concludes that the sequence (^^)neN is bounded in Lfg^{M. x 
R). On the other hand, ^ = (^Ln^^ = ^^n^- Owing to (1.4) and (3.28), the sequence 
(^^)„gN is bounded in Lf^c(R x R). Consequently, up to extraction of another subsequence, 
there exists Vq G Hl^ci^ x R) such that vl^ vo strongly in -^^^^(R x R) and 
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However, a^^ ^ < a ^ >a=< a >h weak-* as n — > +00. Thus, 



dx ar < a>H 



weakly in Lf^^^ x M) as n ^ +00. 



By uniqueness of the hmit, one gets Vq =< a >h Passing to the hmit as n ^ +00 in 
the first equation of (3.22) with L = imphes that Uq is a weak solution of the equation 

~dt " <I^>A 9iuo) =< a >H <I^>A giuo) m V'{R x R). 

From parabolic regularity, the function uo is then a classical solution of the homogenous 
equation 

^<a>H -Q^+ < f^>A g{uo) m K X M, 

such that < 1*0 < 1 and ^ > in R X R. Lastly, JJ^^ Uo{t, x) dtdx^\ from (3.23). On 
the other hand, it follows from the second equation of (3.22) and (3.27) that 

V7eR, uoit+^.x) =UQ{t,x + ^) inMxR, 

where c = lim inf /,_^.o+ = lim„_^.+oo c^„ > 0. In other words, UQ{t,x) — Uq{x + ct), where 
Uq is a classical solution of the equation 

cUo=< a>H Uq+ < fi>Ag{Uo), 0<t/o<linR (3.29) 

that satisfies t/g > in R and 

Standard elliptic estimates imply that Uo converges as s — > ±00 in Cfg^{R) to two constants 
Uq e [0, 1] such that < n >a giU^) = 0, that is giU^) = 0. The monotonicity of Uq 
and the assumption on g imply that Uq = and Uq = 1. In other words, Uo is a usual 
travelling front for the homogenized equation (3.29) with speed c and limiting conditions 
and 1 at infinity. Since the minimal speed for this problem is equal to 2^< a >h< A* >a, 
one concludes that 

lim iiif c*]^ — c> 2-y/^'ar^^<7'^^A- 

4 Monotonicity of the minimal speeds c}^ near the ho- 
mogenization hmit 

This section is devoted to the proof of Theorem 2.3. Before going further in the proof, we 
recall that for each L > 0, the minimal speed is given by the variational formula 

. k{X,L) k{Xl,L) 
Cr = mm — - = — -, 
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where > and k{X,L) is the principal eigenvalue of the eUiptic equation (3.2). Notice 
that /c(A, L) can be defined for all A e R and L > 0. 

Step 1: properties of k{X,L) and definition of k{X,L). The principal eigenfunction 
of (3.2) is L-periodic, positive and unique up to multiplication. Denote 

for all L > 0, A e R and x e R. Each function is l-periodic, positive and it is the 
principal eigenfunction of 

{a(P'^A' + 2LXa(P'^^L + LXa'(f)x,L + L^X''a(f)x,L + L^^i(t>KL = L^k{X, L^x^l, 

associated to the principal eigenvalue L'^k{X,L). But the above problem can be defined for 
all A G M and L G M. That is, for each (A, L) G R^, there exists a unique principal eigenvalue 
k{X,L) and a unique (up to multiplication) principal eigenfunction (f){X,L) of 

{a^'xA' + 2LAa0l_i + LXa'^x,L + L^X''a^x,L + L^I^4>x,l = k{X, L)^x,l. (4.1) 



/' 

Jo 



Furthermore, 0a,l is l-periodic, positive and it can be normalized so that 

' ^lAx)dx = 1 (4.2) 

'0 

for all (A, L) G R^. By uniqueness of the principal eigenelements, it follows that 

VL>0, VAgR, ~k{X,L) = L^k{X,L) 

and 0A,L and 0a,l are equal up to multiplication by positive constants for each L > and 
A G R. ' 

Some useful properties of k{X, L) as L — > O"*" shall now be derived from the study the 
function k. Notice first that, since the coefficients of the left-hand side of (4.1) are analytic 

in (A, L), the function k is analytic, and from the normalization (4.2), the functions 0a,l also 
depend analytically in Hf^M.) on the parameters A and L (see [10, 20]). In particular, the 
function k is analytic in R x (0, -|-oo). Observe also that 

^(A, 0) = and ^a,o = 1 for all A G R. 

Lastly, when A is changed into —A or when L is changed into —L, then the operator in (4.1) 
is changed into its adjoint. But since the principal eigenvalues of the operator and its adjoint 
are identical, it follows that 

V (A, L) G R^ ^(A, L) = k{X, -L) = k{-X, L). 

In particular, it follows that 

r)^k B^ff^'^^k 
V(*,j)eN^ _(A,0) = g^^,^(A,0) = 0. (4.3) 



Therefore, for all A e R, 



k{X,L) = 



~k{X,L) 
L2 



^x0(A,O)as(A,L)^(A,O+). 



But since this hmit is equal to A;(A) = A < a >h -\- < h >a from Step 2 of the proof of 
Theorem 2.1, one then gets that 

1 d'^k — —2 — 

-X— (A,0) = A <a>i^ + </x>A for all A e M. (4.4) 

Z UlJ 

It also follows from (4.3) that 

§(A,L) = ^ X §(A,L) ^ 1 X Jl^ftO) as (A,L) ^ (A,0-). (4.5) 

Prom (4.4) and (4.5), one deduces that 

( A, L) ^ 2 < a >H > as (A, L) ^ (A, 0+) . (4.6) 



9A2 



Similarly, as (A, L) — >• (A,0"'"), 



9A; 9 k{\L) 



d^k 
dXdL 



d^k 



d_ 



' x^(A,L) 



dX 



k{X,L) 



1 d'^k , 



1 



1 ^''^ .X 

X ^(A,0) 



12 dL^ 



(4.7) 



Rerricirk 4.1 As a byproduct of the fact that k and k are even in A, it follows that the 
minimal speed of pulsating fronts propagating from right to left (as in Definition 1.2) is the 
same as that of fronts propagating from left to right. 

Step 2: properties of and A^ in the neighbourhood of L = 0+. Let us first prove that, 
for each fixed L > 0, the positive real number A^ > given in (3.1) is unique. Indeed, if 
there are < Ai < A2 such that 



A;(Ai,L) k{X2,L) 



Ai 



A. 



mm 



k{X,L) 



A>0 A 



then A:(A, L) = c^A for all A G [Ai, A2] since k is convex with respect to A. Then A:(A, L) = c*^X 
for all A G M by analyticity of the map M 9 A 1— k{X,L). But k{0,L) = — > 0, which 
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gives a contradiction. Therefore, for each L > 0, is the unique minimum of the map 

(0,+oo) 3 k{\,L)/\. 

Furthermore, wc claim that L i-^ and L i— > are of class C°° in a right neighbourhood 
of L = 0. Indeed, by definition, A^ satisfies 

F(A* , L) -(A* , L) X A* - k{Xl, L) = 0. (4.8) 

The function (A, L) ^ F(A, L) is of class C°° on M x (0, +00) and |f (A, L) = g|(A, L) x A. 
But 



\* = \J <a>^<li>A > as L ^ 0+ 

from Remark 3.2, and 



dX- 



;{Xl,L) 2 <a>H > as L 



from (4.6). Therefore, from the implicit function theorem, the map L 1— » A^ is of class 
in an interval (0, Lq) for some Lq > 0. As a consequence of formula (3.1), the map L 1— > 
is also of class C°° on (0, Lq). 
For each L e (0, Lq), one has 

by definition of A^ and formula (3.1). But A^ ^ A* > and ||(A2, L) — »■ as L — »■ 0+ from 
(4.7). Thus, 

^ ^ as L ^ 0+. 
On the other hand, it follows from (4.6), (4.7) and (4.8) that 
dX*r. 1 fdk,^^ _ d^k 



Therefore, 



^ — X ^(A*,0) asL^O+, 
12A* 5L4^ ' ^ 

(4.9) 

from (4.7). 

Step 3: calculation of |^(A*,0). In this step, we fix A* = \J^Ca>^l^^i>A- Since the 
functions 0a*,l depend analytically on L G M in Hf^jM.), the expansion 

4>\\L = 1 + + L'^(t>2 + L^(t>z + i-V4 + ■ ■ ■ 
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is valid in Hfg^{M.) in a neighbourhood of L = 0, where 1 = (f)^*^ and 



L=0 



for each i > 1. We now put this expansion into 
and remember that 



and 



MA*,0) = — (A*,0) = — (A*,0) = 



(A*,0) = 2x [{X*f <a>H + <l^>A] =4</x>^ 



from (4.3) and (4.4). Since both 0a*,l and k{X*,L) depend analytically on L, it follows in 
particular that 



' (a(/)'i)' + A*a' = 0, 

(a^'J' + 2\*a(f)[ + \*a'(f)i + {\*Ya + fi = 2 <fi>A, 
{a(j)'J + 2X*a(f)'2 + A*a>2 + {X*fa(j)i + fjL(j)i ^ 2 < h>a 0i, 

(a(f)'J + 2A*a0^ + A*aV3 + (X*yah + //02 = 2 <//>a 02 + x |-^(A*, 0) 



(4.10) 



24 dL^ 

in M. Furthermore, each function 0^ is 1-periodic and, by differentiating the normalization 
condition ||0A*,L||i2(o i) ~ ^ with respect to L at L = 0, it follows especially that 



/ 01 = and / 02 = / 
Jo Jo ^ Jo 



and 

It is then found that, for all x e M, 

0i(a;) = A* X \ -x-\- <a>H / -j- 

V Jo a{l 



dy 



and 

0'2(^) 



<II>A X 



X 



a{x) 



Jo «(y) 



<a>H 



-dy - 



dy - 



<a>H 



<a>H 



dy 



a{x) Jo a{y) 



-dy 



fi{z)dz dy 



fi{y)dy 



+ {yy{x + 



Moreover, it follows from the third equation of (4.10) that, for all x e M, 

a(x)0'3(x) = -2A* ra{y)cf>',{y)dy - A* f d {y)(t>2{y)dy 
Jo Jo 

px px PX 

-{X*f a{y)(j)i{y)dy - ii{y)(pi{y)dy + 2 <fi>A / (t)i{y)dy+ <a>H c, 
Jo Jo Jo 
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where 



= / X fsA* r a(zU'Jz)dz + X* C a' iz)(t)2{z)dz + r a(zUi(z)dz 

Jo [a{y) \ Jo Jo Jo 



+ / fj,{z)(()i{z)dz — 2 < n>A / (l)i{z)dz 



dy 



On the other hand, by integrating the fourth equation of (4.10) over the interval [0, 1], 
one gets that 



{X*,0) = -2<fi>A [ 02 + A*/ a(j>'^ + {yf ! a(j>2+ f fi4>2. (4.11) 
Jo Jo Jo Jo 



Now, put all the previous calculations into (4.11). After a lengthy sequence of integrations 
by parts, it is finally found that 



(A*,0)= / 

^0 



' A{x) 



a{x) 



-dx 



<a>H 



' A{x) 





dx 



where 



A{x) = / i^{y)dy + <ijL>A<a>H / —r^dy - 2<ijl>a x. 
Jo Jo a{y) 

Prom (4.9), it follows that 



dL^ 



7 := 2J <a>H<IJi>A x 



^ A{xf 
a{x) 



dx — <a> 



H 



^ Mx) 

«(2^) 



dx 



as L ^ 0+. 



Cauchy-Schwarz inequality yields 7 > 0. Furthermore, 7 = if and only if A is constant. 
But since ^4(0) = 0, the condition 7 = is equivalent to A'{x) — for all x, which means 
that 

ix{x) , <a>H 



+ 



2 for all X e 



<Ijl>a o,{x) 

In particular, if /i is constant and a is not constant (resp. if a is constant and /i is not 

constant), then this condition is not satisfied, whence limi^o+ -^l^ > in this case. That 
completes the proofs of Theorem 2.3 and Corollary 2.4. □ 

Remark 4.2 In the case when <ij,>a= and ^ 0, then pij^ < for each L > 0, and the 
minimal speed of pulsating traveling fronts is well-defined and it is still positive. Prom 
the arguments developed in this section and in the previous one, one can check that, in this 
case. 







+ 



where 



0+, 



dX*, 



dL V <a>H 

^ A{x)^ 



> and 



a{x) 



dx — <a> 



H 



del _ 
dL 

^ A{x) 
a{x) 



2^/3 <a>H > as L ^ 0+, 



dx] > 
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and A{x) = / fi{y)dy. Therefore, the speeds are increasing in a right neighbourhood 
Jo 

of L = but, in this case, the variation is of the first order. Notice that the formula 
hmi^o+ ~aL ~ 2-\//3 <a>H is coherent with the numerical calculations done by Kinezaki, 
Kawasaki and Shigesada in [21] (see Figure 3b with < // >a— 0, that is A — under the 
notations of [21]). 



5 Proof of Theorem 2.6 

As in the proofs of the previous theorems, we use the following formula for the minimal 
speed: 

c* = mm — ^ = — (5.1) 
^ A>o A A* ^ ' 

where kz{\) is defined as the unique real number such that there exists a positive Lo-periodic 
function -0 satisfying: 

< + 2A V'' + A'V + l^z{x)^l: = A;,(A)^ in (0, Lq). (5.2) 
Setting (p{x) = e^^ip{x), the above equation and periodicity conditions become equivalent 

to: 

(p" + Hz{x)^p = in (0, Lo), 

V9(Lo) = e^^V(O), (5.3) 
(^'(Lo) = e^^0(^'(0), 

which therefore admits, for every positive A, a unique solution ((/?, kz{\)) with (/? > satisfying 
the normalisation condition (/^(O) = 1. 

Let A > be fixed. System (5.3), together with the normalization condition (p{Q) = 1, is 
equivalent to: 

ip" = {kz{X) -m)<f on [0,1/2), 
ip" = kz{X)p on [1/2,1/2 + z), 

^" = {kz{X) -m)ip on [1/2 + z, I + z), (5.4) 
ip" = kz{X)ip on [I + z,Lo), 
(<p{0) = 1, (/.(Lo) = e^^V(0), (^'(Lo) = e^^O(/.'(0). 

For each z e [0, Lq — /], let A* be defined by the formula (5.1). We have the following 
lemma: 

Lemma 5.1 Assume that I > 3Lq/4. Then, for all z e [0, Lq — I], we have kz{Xl) > m. 

Proof of Lemma 5.1. Let us divide equation (5.2) by '4> and integrate by parts over [0, Lq]. 
Using the Lo-periodicity of ^0, we obtain: 

»Lo u/,/|2 /-Lo 



/ ^ + LoA2+ / fXz{x)dx^ Lokz{X). 
Jo W Jo 
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Thus, 

1 f^" I 
A:^(A)>A^ + — / = + m— . (5.5) 

Lo Jo Lo 

From (5.1) and (5.5) we get: 

Thus, (A*)^ - 2Xly/m + ml/Lo < 0, which imphes that 

A* > -\A^ ~ yjm — ml/Lo. 

Using (5.5), we finally get 

k,{Xl) > 2m(l - y/l - //Lo) > m, 

as soon as I > SLq/A. □ 
We now turn to the proof of Theorem 2.6 and we assume that / G (3Lo/4, Lq). Using 
the fact that ip G C^(M), a straightforward but lengthy computation shows that, whenever 
kz{X) > m, system (5.4) is equivalent to 

F{z,X,kz{X)) _ 
G{z,k,{X)) 

where F and G are two functions, defined respectively in [0, Lq — I] x (0, +oo) x [m, +oo) 
and [0, Lq — I] X [m, +oo) by: 

F{z, A, s) = 4(2s — m)^/s^/s~^^msmh{l^/s~^^rn) sinh(Q;y^) 

cosh.{(3^/s){l — cosh.{l^/s — m)) , . 

+8(s^ — ms)[cosh.{ly/s — ni) cosh(ay^) — cosh(ALo)] 
cosh(Q;-^)(cosh(/-\/s — m) — 1), 



and 



G{z^s) — m ^/s cosh{l y/s~^^rn) [4 sinh(a; y^)(s/m — 1) 

+ (sinh(Q;y^) — smh{P^/s)) (l — 1/ cosh(/Vs — m))] 



+m^/s^^msmh.{l^/s^^rn) cosh(Q;y^) 



4s _ , cosh(/3^) 
m cosh(a\/s) 



(5.7) 



with a := Lq — I and (5 :— Lq — I — 2z. 

Each factor in the expression (5.7) is positive, as soon as s > m, for z G [0, Lq — I]. Thus, 
whenever ^^(A) > m, system (5.4) is equivalent to the simpler equation 

F{z,X,k,{X))^0. (5.8) 

Furthermore, from Krein-Rutman theory, since the eigenfunction tjj in (5.2) is positive, kz{X) 
is the largest real eigenvalue of the operator t/j t/j" + 2A ip' + A^V + IJ'z{x)'^. This result, 
implies that, for each z G [0, Lq — I], and each A > 0, kz{X) is the largest real root of equation 
(5.8), as soon as kz{X) > m. 
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Prom equation (5.6), we easily see that 

lim F{z,X,s) = +00, (5.9) 

s— ++00 

for all z e [0, Lq — I] and A > 0. Moreover, differentiating (5.6) with respect to z, we obtain 

— (2;, A, s) — 2m^-\/ssinh(-\/s(Lo — / — 2z)) [cosh.{l^/s — m) — l] . 

Thus, for all s > m, and A > 0, 

dF 

— {z, X,s)>0 for z e [0, (Lo - /)/2), (5.10) 
oz 

and 

dF 

— {z,X,s)<Oiorze{{Lo-l)/2,Lo-l]. 

Now, take Zi < Z2 in [0, {Lq — I)/ 2], and assume that c*^ < c*^. It follows from formula 
(5.1) that kz2{X) > cl^X, for all A > 0. In particular. 

Prom Lemma 5.1, we know that A;2^(A*J > m. Thus, (5.11) implies A;22(A*J > m. 
From the above discussion, /c^2(A*J is therefore the largest real root of the equation 
F{z2, A*^, A;^2(A*J) = 0, and, similarly, A'"^i(A*J is the largest real root of F{zi, A*^, /c2i(A*J) = 
0. Using (5.9) and (5.10), and since < zi < Z2 < {Lq - 1)/2, we obtain kz2{XlJ < fczi(A*J, 
which contradicts (5.11). Therefore, c* is a decreasing function of z in [0, (Lq — /)/2]. Similar 
arguments imply that c* is an increasing function of z in [(Lq — l)/2, Lq — I]. This concludes 
the proof of Theorem 2.6. □ 
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